Characterization of ferrimagnetic Heisenberg chains according to the constituent spins 
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The low-energy structure and the thermodynamic properties of ferrimagnetic Heisenberg chains 
of alternating spins 5* and s are investigated by the use of numerical tools as well as the spin- 
wave theory. The elementary excitations are calculated through an efficient quantum Monte Carlo 
technique featuring imaginary-time correlation functions and are characterized in terms of interact- 
ing spin waves. The thermal behavior is analyzed with particular emphasis on its ferromagnetic 
and antiferromagnetic dual aspect. The extensive numerical and analytic calculations lead to the 
classification of the one-dimensional ferrimagnetic behavior according to the constituent spins: the 
ferromagnetic {S > 2s), antiferromagnetic (5* < 2s), and balanced (5* = 2s) ferrimagnetism. 

PACS numbers; 75.10.,]m, 65.50.+m, 75.40Mg, 75.30Ds 



I. INTRODUCTION 

One-dimensional quantum ferrimagnets are one of the 
hot topics and recent progress [p}^15[ in the theoretical 
understanding of them deserves special mention. Such 
a vigorous argument a great deal originates in the pio- 
neering efforts to synthesize binuclear magnetic materials 
including one-dimensional systems. The first ferrimag- 
netic chain compound MnCu(dto)2(H20)3-4.5H20 
(dto = dithiooxalato = S2C2O2), was synthesized by 
Gleizes and Verdaguer and stimulated the public interest 
in this potential subject. The following example of an 
ordered bimetallic chain MnCu(pba)(H20)3-2H20 
(pba = 1, 3-propylenebis(oxamato) = C7H6N2O6), ex- 
hibiting more pronounced one dimensionality, further ac- 
tivated the physical [Il8|jl9[ , as well as chemical in- 
vestigations. There also appeared an idea that the 
alternating magnetic centers do not need to be metal ions 
but may be organic radicals. 

A practical model of a ferrimagnetic chain is two kinds 
of spins S and s {S > s) alternating on a ring with anti- 
ferromagnetic exchange coupling between nearest neigh- 
bors, as described by the Hamiltonian, 



N 



(1) 



Here N is the number of unit cells and we set the unit- 
cell length equal to 2a in the following. The simplest 
case, {S, s) = (1, i), has so far been discussed fairly well. 
There lie both ferromagnetic and antiferromagnetic long- 
range orders in the ground state The ground state, 
which is a multiplet of spin {S — sYN, shows elementary 
excitations of two distinct types M. The excitations of 



ferromagnetic aspect, reducing the ground-state magne- 
tization, form a gapless dispersion relation, whereas those 
of antiferromagnetic aspect, enhancing the ground-state 
magnetization, are gapped from the ground state. As 
a result of the low-energy structure of dual aspect, the 
specific heat shows a Shottky-like peak in spite of the 
ferromagnetic low-temperature behavior and the mag- 
netic susceptibility times temperature exhibits a round 
minimum |^,|l^,|l8|. When the exchange coupling turns 
anisotropic, the dispersion of the ferromagnetic excita- 
tions is no more quadratic [Q and the plateau in the 
ground-state magnetization curve due to the gapped 
antiferromagnetic excitations vanishes via the Kosterlitz- 
Thouless transition |]l5| . 

The quantum behavior of the model with higher spins 
is not yet so clear as that in the {S,s) = (1, |) case. 
Though Drillon et al. 18 1 made the first attempt to re- 
veal the general behavior of the model, they fixed the 
smaller spin to ^ in their argument with particular em- 
phasis on the problem of the crystal engineering of a 
molecule-based fcrromagnet— the assembly of the highly 
magnetic molecular entities within the crystal lattice in 
a ferromagnetic fashion. There also exists an extensive 
numerical study but the leading attention was not 
necessarily directed to the consequences of the variation 
of the constituent spins. So far the generic behavior, 
rather than individual features, of ferrimagnetic mixed- 
spin chains has been accentuated or predicted, but there 
are few attempts to characterize or classify the typical 
one-dimensional ferrimagnetic behavior as a function of 
{S,s). In such circumstances, we aim in this article at 
elucidating which property of the model (M is universal 
and which one, if any, is variable with {S, s) . 

It is true that numerical tools are quite useful in this 
context, but they are not almighty. Although the low- 
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temperature ferromagnetic behavior is quite interesting 
from the experimental point of view, it is hardly feasi- 
ble to numerically take grand-canonical averages at low 
enough temperatures. It is also unfortunate that with 
the increase of (5, s), the available information is more 
and more reduced in both quality and quantity. Then 
we have an idea of describing the model in terms of the 
spin-wave theory. The conventional spin-wave treatment 
of low-dimensional magnets may discourage us. The Hal- 
dane conjecture revealed a qualitative breakdown of 
the spin-wave description of one-dimensional Heisenberg 
antifcrromagncts. Neither quantum corrections nor 
additional constraints to spin fluctuations end up 
with an overall scenario for the low-energy physics ap- 
plicable to general spins. However, Ivanov [Q has re- 
cently reported that the spin-wave description of one- 
dimensional Heisenberg ferrimagnets is quite successful. 
Though his calculations were restricted to the ground- 
state energy and magnetization, the highly accurate es- 
timates obtained there are surprising enough, consid- 
ering the diverging ground-state magnetization in the 
one-dimensional antiferromagnetic spin- wave theory. For 
antiferromagnets, quantum fluctuations of domain-wall 
type, connecting the two degenerate Neel states, are im- 
portant, whereas for ferrimagnets, domain-wall excita- 
tions lead to magnetization fluctuations and are thus of 
less significance. Therefore it is likely that the spin-wave 
approach is highly efficient for ferrimagnets. In an at- 
tempt to demonstrate such an idea, we try to describe 
not only the low-energy structure but also the thermal 
behavior of the Heisenberg ferrimagnetic spin chains (|l|) 
within the framework of the spin-wave theory. Extensive 
numerical calculations, supplemented by the spin-wave 
analysis, fully reveal the one-dimensional ferrimagnetic 
behavior as a function of the constituent spins. 



II. ELEMENTARY EXCITATIONS 

In order to investigate the low-energy structure, we em- 
ploy a new quantum Monte Carlo technique pst as well as 
the conventional Lanczos diagonalization algorithm. The 
idea is in a word expressed as extracting the low-lying 
excitations from imaginary-time quantum Monte Carlo 
data at a low enough temperature. The imaginary-time 
correlation function S{q, r) is generally defined as 



5(g,T)~^|(l;fco|5||/;fco + g)| 



'T[Ei(ko+q)~Ei{ko)] 
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where Oq — N^^^-^^Oje'^^"-'^^ is the Fourier trans- 
form of an arbitrary local operator Oj, which may be 
an effective combination of the spins S and s, and the 
thermal average at a given temperature /3~^ = k^T, 
{A) = Tr[e"'^'^A]/Tr[e"'^'^], is taken in a certain sub- 
space. S{q,T) can be represented in terms of the eigen- 
vectors and eigenvalues of the Hamiltonian, |/; fc) {I — 
1,2,- • •) and Ei{k) {Ei{k) < E2{k) < ■■■), and behaves 
like 



(3) 



at a sufficiently low temperature, where kg is the mo- 
mentum at which the lowest-energy state in the sub- 
space is located. Therefore we can reasonably approx- 
imate i?i(fco + q) — Ei{ko) by the slope — i91n[S'(g, t)]/9t 
in the large-r region. When we take interest in the 
lower edge of the excitation spectrum, such a treatment 
is rather straightforward. The elementary excitations of 
the Haldane antiferromagnets were indeed revealed thus 
p6[ . Here, due to the two kinds of spins in a chain and 
the resultant dual aspect of the low-energy structure, the 
relevant subspace and operator Oj are not uniquely de- 
fined. Since the total magnetization, M = J^ji'^j + ^j)^ 
is a conserved quantity in the present system, we con- 
sider calculating S{q,T) independently in each subspace 
with a given M |27|. The elementary excitation energies 
for the ferromagnetic branch are obtained from a single 
calculation, S{q, r) in the subspace of Af — 0, while those 
for the antiferromagnetic branch from a couple of calcu- 
lations, S{q, t) and the lowest energy in the subspace of 



M = {S- s)N 
found that O 



h 1. We have taken 



± Sj for Oj 



and 



J — ^j — Sj extracts the eigenvalues of the 
bonding (lower-energy) states in both subspaces. The 
choice of the scattering matrices is a profound problem 
in itself and is fully discussed elsewhere P,|2^. 

Although the chain length we can reach with the exact- 
diagonalization method is strongly limited, the diago- 
nalization results are still helpful in the present system 
whose correlation length is generally so small as to be 
comparable to the unit-cell length (see Fig. || bellow). 
Actually, the ground-state energies for = 10 coincide 
with their thermodynamic-limit values within the first 
several digits. The Lanczos algorithm gives the most pre- 
cise estimate for the ground-state energy and the antifer- 
romagnetic excitation gap, whereas the quantum Monte 
Carlo technique is necessary for the evaluation of the cur- 
vature of the small-momentum dispersion. 

On the other hand, we consider a spin-wave descrip- 
tion of the elementary excitations as well. We introduce 
the bosonic operators for the spin deviation in each sub- 
lattice via 



Sj' = y2S — a]aj , Sj = S — aja^ , 



(4) 



where we regard S and s as quantities of the same order. 
The Hamiltonian (|l|) is expressed in terms of the bosonic 
operators as 



H^E,uss+no + Hi + 0{S-^). 



(5) 



where E'ciass = ~2sSJN is the classical ground-state en- 
ergy, and Ho and Hi are the one-body and two-body 
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5uof = 2(5* + s) A 



rf + ri + (^^s + ./7/s) , (11) 

sin^(afc) 



[iOk ±{S~ s)] , (12) 



with 



k 



i^Ss cos^(afc) 



(13) 



while the irrelevant one-body terms 



^irrcl J~ 
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aE^k/^^ + 44), (14) 



and the residual two-body interactions Tiquart are both 
neglected so as to keep the ferromagnetic branch gapless. 
The resultant Hamiltonian is compactly represented as 



with 



= uj, — Suj, 



Eg; = -Eclass + Eq 
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(15) 
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Now we compare all the calculations in Fig. |^. The 
diagonalization results fully demonstrate the steady ap- 
plicability and good precision of our Monte Carlo treat- 
ment. The spin-wave approach generally gives a good 
description of the low-energy structure. Even the free 
spin waves allow us to have a qualitative view of the 
elementary excitations. The relatively poor description 
of the antiferromagnetic branch by the free spin waves 
implies that the quantum effect is more relevant in the 
antiferromagnetic branch. Here we may be reminded of 
the spin-wave treatment of mono-spin Heisenberg mag- 
nets. For the ferromagnetic chains, the spin-wave disper- 
sions are nothing but the exact picture of the elemen- 
tary excitations, while for the antiferromagnetic chains, 
those are generally no more than a qualitative view. The 
point is that in the present system the spin-wave picture 
is efficient for both elementary excitation branches and 
the spin-wave series potentially lead to an accurate de- 
scription. More specifically, the divergence of the boson 
numbers, which plagues the antiferromagnetic spin-wave 
treatment in one dimension, does not occur in the present 
system. This viewpoint is further discussed in the final 
section. We conclude this section by listing in Table || 
the spin- wave estimates of a few interesting quantities in 
comparison with the numerical solutions, where we define 
the curvature v as ojf^^Q — v{2ak)^. 



III. THERMODYNAMIC PROPERTIES 

The dual structure of the excitations leads to unique 
thermal behavior. In order to complement quantum 
Monte Carlo thermal calculations, especially at low tem- 
peratures, we consider describing the thermodynamics 
in terms of the spin-wave theory. Introducing the addi- 
tional constraint of the total magnetization being zero 
into the conventional spin-wave theory, Takahashi p9| ] 
succeeded in obtaining an excellent description of the 
low-temperature thermal behavior of one-dimensional 
Heisenberg ferromagnets. The present authors have re- 
cently applied the idea to the spin-(l, i) ferrimagnetic 
Heisenberg chain |jl^. Here we develop the method for 
general spin cases and make a detailed analysis of its 
vahdity as a function of {S, s). The core idea of the so- 
called modified spin- wave theory |30 3^ can be summa- 
rized as controlling the boson numbers by imposing a 
certain constraint on the magnetization. From this point 
of view, the zero-magnetization constraint, which is quite 
reasonable for isotropic magnets, plays a relevant role in 
ferromagnets. The resultant low-temperature expansions 
pof of the specific heat and susceptibility of the spin-s 
Heisenberg ferromagnetic chain. 
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with t = k^T/J and Riemann's zeta function C{z), in- 
deed coincide with the thermodynamic Bethe-ansatz cal- 
culations for s = 5 within the leading few terms. 

In the present system, the zero-magnetization con- 
straint is explicitly represented as 

E(5| + 4) = iS-s)N-J2Y. ^^k" = ' (20) 



where 



1+) with Pfc(n ,n'^) be- 



ing the probability of n~ ferromagnetic and antiferro- 
magnetic spin waves appearing in the fc-momentum state. 
Equation (2^) straightforwardly proposes that the ther- 
mal spin deviation should cancel the Neel-state magneti- 
zation. By minimizing the free energy 

k 

+kBTj2 E Pkin-,n+)\nPk{n-,n+), (21) 

k n^,n+ 

with respect to Pfc(n~,n+) at each k under the 
condition ( [20| ) as well as the trivial constraints 
J2n- n+ Pk{n^ ,n+) = 1, we obtain 
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t-i + 0(t-5), (23) 



where i = ksT/Jj with 7 = 1 ^ ri{S + s)/ Ss - r2/\^. 
The specific heat C has been obtained by differentiating 
the free energy F, whereas the susceptibihty x by calcu- 
lating ((M^) — (M)^)/3r, where we have set the g factors 
of the spins S and s both equal to g, because the differ- 
ence between them amounts to at most several percent 
of themselves in practice Q . 
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FIG. 2. Schematic representations of the M — {S — s)N 
ground states of spin-(S, s) ferrimagnetic chains of A'' elemen- 
tary cells in the decoupled-dimer limit (a) , the AKLT ground 
states of spin-(2s) antiferromagnetic chains of 2N spins (b), 
the M = 2N{S — s) ground states of spin-(S' — s) ferromag- 
netic chains of 2A'' spins (c). The arrow (the bullet symbol) 
denotes a spin 1/2 with its fixed (unfixed) projection value. 
The solid segment is a singlet pair. The circle represents an 
operation of constructing spins S and s by symmetrizing the 
spin 1/2's inside. The relation (a) « [(b) -I- (c)]/2 may be 
expected. 



The analytic expressions (|2^) and (|2^) give us a bird's- 
eye view of the one-dimensional ferrimagnetic behavior. 
The spin-(5', s) ferrimagnet turns into the spin-s antifer- 
romagnet in the limit of S' — > s, whereas it looks like 
the spin-5' ferromagnet in the limit of S/s — > 00. In 
this sense, the subtraction S — s may be regarded as 
the ferromagnetic contribution, while the residual spin 
amplitude 2s as the antiferromagnetic one. No ferro- 



magnetic aspect at S/s = 1, while a hundred percent 
ferromagnetic aspect for S/s — > 00. Another consider- 
ation also leads us to such a picture. Since the pertur- 
vation from the decoupled dimers jsj qualitatively well 
describes the low-lying excitations of the system, we pro- 
pose in Fig. ^ an idea of decomposing ferrimagnets 
into ferromagnets and antiferromagnets, where we let the 
decoupled dimers and the Affleck-Kennedy-Lieb-Tasaki 
valence-bond-solid states 34 1 symbolize ferrimagnets and 
integer-spin gapped antiferromagnets, respectively. Now 
we expect spin- (5,5) ferrimagnets to behave like combi- 
nations of spin-(2s) antiferromagnets and spin-(5 — s) 
ferromagnets. Since the antiferromagnetic excitations 
of the ferrimagnetic ground state are gapped, the low- 
temperature behavior of ferrimagnets should be only of 
ferromagnetic aspect. At low temperatures there is in- 
deed no effective contribution from the spin- (2s) antifer- 
romagnetic chain with an excitation gap [ p2[ immediately 
above the ground state. In this context, we are surprised 
but pleased to find that provided S = 2s, the expres- 
sions ( p^ ) and ( p3| ) coincide with those for ferromagnets, 
( p^ ) and ( p^ , except for the quantum renormalizing fac- 
tor 7. The low-temperature thermodynamics should be 
dominated by the small-momentum ferromagnetic excita- 
tions. Within the linearized spin- wave theory, the small- 
momentum dispersions of Heisenberg ferrimagnets and 
ferromagnets are characterized by the curvatures 
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J, 



2(5' -s) 



(24) 
(25) 



respectively, where we find that they coincide with each 
other only when S = 2s. The criterion, S — 2s, is fur- 
ther convincing when we consider the high-temperature 
behavior. The paramagnetic behavior of the spin-(S', s) 
ferrimagnet is given as 



-fciT 



NksT 



= -ln[(25 + l)(2s + l)] , 



(26) 



\% 2 = o iS{S + 1) + s(s + 1)] , (27) 

whereas those of the spin-(5 — s) ferromagnet and the 
spin-(2s) antiferromagnet are as follows: 

j^{S — s)-icTvo ^(2s)-antifcrro 

NksT 

= -ln[(25'- 2s-|- l)(4s-|- 1)] , (28) 

^gy(-Y(5'-s)-fcrro _|_ ^(2s)-antifcrro^ 

= i[(5-s)(S'-s + l) + 2s(2s + l)] . (29) 
o 

These asymptotic values agree with each other only when 
S = 2s. This is simply the consequence of the degrees 
of freedom. In ferrimagnets of S* > 2s, the ferromagnetic 
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TABLE I. Comparison between the linear-spin- wave (LSW), interacting-spin-wave (ISW), and numerical estimates of the 
ground-state energy Eg, the antiferromagnetic excitation gap A, and the curvature of the ferromagnetic dispersion, v. The 
most accurate numerical estimates of Eg and A are obtained by the exact diagonalization (Exact), whereas those of v by the 
quantum Monte Carlo (QMC). 



iS,s) 




Eg/NJ 






A/J 






V 




LSW 


ISW 


Exact 


LSW 


ISW 


Exact 


LSW 


ISW 


QMC 




-1.4365 


-1.4608 


-1.4541(1) 


1 


1.6756 


1.759(1) 


1 

2 


0.3804 


0.37(1) 


V 2 ' 2 / 
(1,1) 

(2,1) 


-1.9580 


-1.9698 


-1.9672(1) 


2 


2.8025 


2.842(1) 


3 


0.3390 


0.31(1) 


-3.8281 
-4.8729 


-3.8676 
-4.8973 


-3.861(1) 
-4.893(1) 


1 
2 


1.5214 
2.6756 


1.615(5) 
2.730(5) 


2 
1 


1.1319 
0.8804 


0.90(3) 
0.76(2) 
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